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Abstract 

It was shown that an infinite sequence of improving non-increasing upper bounds to the 
ground state energy (GSE) of a slow-moving piezoeletric polaron can be devised. 
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1. The piezoelectric polaron model 


A quantized polaron model for the case of an electron interacting with acoustic phonons 
through piezoelectric deformation potential was introduced by A.R. Hutson [TJ, its Hamiltonian 
structure being similar to the Frohlich optical polaron model introduced by H. Frohlich [3] 
earlier: 


(1) H = f- + £ fcki-^k + E ft fee-" + 6 k e' 

k k 

where Wk = -sk is the frequency of the acoustical phonons with s being the velocity of sound, 


where V is the volume of the crystal, and 


m 


1 2 <4fc> 

ry = —p - ; - 

2 eCsh 

is the dimensionless coupling constant where is an average of the piezoelectric tensor [l], 

e is the dielectric constant and C is an average elastic constant. The operators p and r are the 
electron momentum and position coordinate quantum operators, 

\Pi,rj\ = 

and the Bose operators & k , b k, 


[^ki ^k'] — ^kk'j [^ki ^k'] — 0; 

create and annihilate phonons of one effective acoustic mode of energy hw k and wave vector k. 

In what follows the energies will be expressed in units of 2 ms 2 , the length in units of h/2ms 
and the phonon wave vectors in units of 2ms/h so that all variables are dimensionless. In this 
units the model CD) reads as 


( 2 ) 


with 


H = P 2 + Y, k K h * + E V * ( 6 k e "' kr + ^e ikr ) , 

k k 


where V is dimensionless volume. Ultimately, the sum over the phonon vectors is to be 
replaced by the integral U/(27r) 3 f dk with a finite cutoff at k o, the boundary of the hrst Brillouin 
zone in the phonon wave vector space, which is introduced to account for the discreteness of the 
crystal lattice with k o ~ 1/a, where a is the lattice constant. 





2. Piezoelectric polaron GSE 


As usual, one can introduce the polaron total momentum 

P = P + E k6 k 6k 

k 

being a constant of the motion and commuting with the Hamiltonian dU). Hence, the Hamilton¬ 
ian can be transformed to the representation in which P becomes a ”c”-number P, the value of 
the total polaron momentum, 

H —»• H, H = S~ l HS, S = exp(—i kr6^b k ), 

k 

(3) H = (P-Yl k6 k 6 k) 2 + X] k K b * + E Vk ( 6 k + bk) , 

k k k 

and the Hamiltonian ([3]) does not contain the electron coordinates anymore. One more subse¬ 
quent transformation 


H —t 'H(f), 

results in the Hamiltonian 


H(f) = U~ 1 HU, t/ = exp{^/ k (6+-6 k )}, 

k 


(4) H(f) = H 0 (f) + H 1 (f), 

where 


Ho(f) = P 2 + E k K b * + (E k6 k ft k) 2 - a', 

k k 

Hi (/) = EP + fc2 )/k + Vk\K + 6 k ) + 2^(k • m)/ k / m 6+6 m + 

k km 

+B k - m )/k/m(6 k 6+ + bkbm) + 2 E< k - m)/ k (6+6 m 6 k + 6+6+ 6 m ) - 

km km 

-2^(P-k)(6+ + / k )(6 k + / k ) + 

k 

(5) +2 E^ k ‘ m )/m b k b k + 2 E(k • m )/m(^k + 6 k) + E( k • m )/m/k- 

km km km 

and 


= 2 E ^/k + E( fc + fe2 )/k + (E /k k ) 2 - 

k k k 
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3. On improvable upper bounds to various polaron models GSE in general 


Our objective here is to demonstrate that an infinite, in principle, sequence of improvable 
upper bounds to the ground state energy E(a, P, ko) of the Hamiltonian (JH) can be constructed 
in a regular way by means of a variational method outlined in [2]. Similar variational approach 
was formulated later in [5]. This energy corresponds to the lowest energy of the slow-moving 
polaron for a given value of the total polaron momentum P. Of course, some reservations 
are to be made here regarding the fact that, unlike in the case of optical polaron, there is no 
energy separation between the polaron ground state and the excited states with the same total 
momentum. Therefore, what is considered here is actually a zero temperature polaron behavior. 
Then, expansion of the function E g (a, P, ko) in powers of P of the kind 

p 2 

Eg(a, P, ko) = E g (a, 0, ko) + —-b 0(P 4 ), 

2m eff 

where E g (a, 0, ko) is the GSE of the polaron at rest, may provide us with some approximate 

value of the polaron effective mass m e /y in spatially isotropic as well as anisotropic case, in 

which general case the so-called inverse effective mass tensor I —-— 1 is to be introduced as 

V r l ef f / ij 

( 1 \ _ d 2 E(a,P,fc 0 ) 

dPidPj p=0 ' 

instead of the scalar effective mass parameter m e ff. 

As was shown in [4] , for a quantum system Hamiltonian H and a trial state | ijj), such that 

(Y#) = 1, 


Eg < min(a[ n) ,...,<4 n) ) < (ip\H\ 


where the real numbers (a^ n \ ...,an L> ) are the roots of the n-th order polynomial equation 


, 0 ) 


Pn(x) = Y J X i X n ~ i = 0. 

i=0 

Here, the coefficient Xq = 1 and all the other coefficients X^. 1 < i < n satisfy the system of n 
linear equations 


with 


MX + Y = 0, 


jVl2n—i i Xiij 'If 2// (i - j) ■ Yj 1,2,.. . 77 ., 

and 


Mm = 

It is assumed that all moments M m are finite. It was proved that the following inequality holds 
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So, at the first order 


/ (^+1) (n+l), ^ . / (ra) (77) \ 

mm(ai \...,a ' n+1 ) < mm(a\ ..., E n >). 

(6) E g <a{\ a {3) = (ip\H\ip), 

and at the second order 


( 7 ) 


Eg < min(aS 2 ) , 4 2) ) = (^\H\^) + ^ 7 - 


a i ] = (VI##) + 


= (VI##) + 77T^ + 


K-, 

2Ko 


k 3 


2Ko 


where K 2 and K 3 are the central moments 


K 3 
2K-2 


+ K 2 


1/2 


K 3 

2K, 


#3 

2 Ko 


11/2 


+ I<2 


1 1/2 


+ I<2 


k 2 = #i (h - #i#i#) 2 i#, k 3 = mn - 

It is seen that the second order bound m would lie below the first order bound (JUJ) for most 
physically relevant quantum models and most reasonable and meaningful choices of the trial 
state \ip). 


4. Improvable upper bounds for piezoelectric polaron at rest 

For P = 0, the function /k is spherically symmetric, and the piezoelectric polaron model- 
Hamiltonian ((U) reduces to 


u(f) = J2 kb i b * + E k6 k 6k ) 2 -«' + 

k k 

+ + k 2 )fk + Vk](b^ + &k) + 2^^(k - m)/fc/ m 6^6 m + 

k km 

+ ' m )/fc/"*( b k b m + Mm) + 2 ^(k • m)/ fc (6+ b m b k + &£&+6 m )- 

km km 

As to be seen later on, it would be convenient to choose phonon vacuum state |0) as a trial state 
\ip) for 'H(f), so that inequality 

E g {a, 0, ko) < (0\H(f)\0) = 2^ V k f k + + k 2 )f 2 

k k 

holds, the right-hand side of which is minimized by 

fk = ~V k /(k + k 2 ), 
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and, eventually, 


( 8 ) 


2 a 

E g (a,0,k 0 ) < E w (a,0,k 0 ) = ln[l + k 0 \. 

7T 


The upper bound (J5J) steins from the conventional variational principle in quantum mechanics 
and is valid for arbitrary value of a. To derive a sequence of better non-increasing upper bounds 
one needs to calculate moments (0\'H m (f)\0) for sufficiently large integer exponents m by means 
of the Wick theorem. The result of the calculation for any such moment can be expressed in 
terms of the products of integrals of rational functions 


( 9 ) 


rko 

Jo 


k v dk 


p,q — non-negative integers, 


o (k + k 2 )'?’ 
which can be evaluated analytically. 

Thus, at the second order variational approximation J7J) 


2 a K-> 

(10) E g (a, 0, ko) < E var = ln[l + k 0 ] + —- 

7r Ah-2 

where 


Kj 

2K 2 


+ K 2 


1/2 


( 11 ) 


K-2 = ^Fi(ko), 


( 12 ) 


16a 2 64a 3 o 

K 3 = —(F^F^ko) + Fi(fco)F 3 (A:o)) + —^F?(k 0 ), 


3ir 2 


97T 3 


(13) 

(14) 


Fi(k 0 ) = 


ln(l + fcg) + 


l + ^o 


-1 


F 2 (k 0 ) = 


kn 


- k 0 + ln(l + ko) 


(15) 


Fo{ko) = 


k 2 1 

—2ko + -F + 3 ln(l + ko) + -—— 

2 1 + k 0 


- 1 


This upper bound is worth comparing with the lower bounds to the piezoelectric polaron GSE 
obtained in [6] for the case of small 


(16) 

and large 


Flbs ~ -(2a/7r)ln(k 0 + 1), a < 1, 


( 17 ) 


Flbl ~ --a 2 — (4a/7r) ln(fe 0 /a), l<aCio 

O 
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coupling constant respectively. The bound Elbs coincides, actually, with the perturbation 
theory result [3- It was assumed in [6] that ko ~ 150. So, the bounds (HD, (USD and m are 
plotted in Figs. 1-3 just for this value of fco- 


5. Improvable upper bounds to the GSE of the slow-moving piezoelectric 

POLARON 


In general case P / 0 


E g (a, P, ko) < <0|P(/)|0> =P 2 + 2 E V k f k + + k 2 )f 2 - 2 £(P • k )/ 2 + (£, /k k) 2 , 

k k k k 

with the right-hand side to be minimized by 

h = -V k /[k-2k-P(l-r,) + k 2 }, 

where t] is defined self-consistently by the equation 

r/P = E /kk = E V kVl k ~ 2k ' ^ ~ ? ?) + 

k k 

or, alternatively, by 


r/P 2 = ifk • P /[k - 2k • P(1 -n) + k 2 } 2 . 

k 

The resulting upper bound is 


(18) 


E g (a,P,k 0 )<P 2 (l-r,) 2 ~J2 V k 

k 


k + k 2 -4k-P(l -??) 

[k — 2k • P(1 — rj) + A; 2 ] 2 ' 


Another choice 


(19) / k = -[Vfc + 2r?k • P]/[As - 2k • P + fc 2 ], 

eliminating all terms linear in Bose operators b^. in Q, is possible too, with the corresponding 
self-consistency equation for 7] 

ri p2 = E f* k = E k ' p \ y k + 2r?k • P } 2 /[k - 2k • P + fc 2 ] 2 , 

k k 

which can be solved analytically. At the same time, the simplest choice 

fk = ~V k /(k + k 2 ) 

seems to be preferable, because in this case the complexity of the analytical calculation of 

arbitrary moments (0|P m (/)|0) does not exceed the one for the case P = 0, i.e. no employment 
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of any integrations over wave vectors more complicated and laborious than the integrals of the 
type ([9|) is necessary. 


6. Summary 

It was shown that the GSE function E g (a, P,&o) of the slow-moving piezoelectric polaron 
can be approximated by infinite sequence of non-increasing upper bounds for arbitrary values of 
the coupling constant a , polaron total momentum P and cut-off wave vector ktj. The proposed 
algorithm for the construction of these bounds is well-suited for implementation by means of 
contemporary techniques for parallel computing due to its overwhelming reliance on the Wick 
theorem. 
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Figure 1. Upper bounds: E var , 
solid line; Elbs, dashed line; 
Elbl-, dash-dotted line; fco = 
150. 
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